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Resonant excitations of confined systems have aroused much attention because of their potential 
appUcation in future microwave devices and spintronics. Under resonant excitations, the motion of 
topo-logical objects exhibits circular, elliptical or even stadium-like dynamics. However, more complex 
non-linear resonant excitations of topological objects in confined systems have seldom been reported and 
the associated physical mechanism remains unclear. Here, we present an observation of flower-like resonant 
excitations for coupled skyrmions in Co/Ru/Co nanodisks activated by a single-frequency microwave 
magnetic field by means of numerical simulation. We find that flower-like dynamics of coupled skyrmions is 
always accompanied by an excitation of an eigenfrequency near 1.15 GHz, which is strongly associated with 
the large non-local deformation of the topological density distribution of coupled skyrmions. These results 
distinguish a skyrmion from other topological objects in dynamics and will be instrumental to the 
manipulation of skyrmions for applications. 

Since proposed by Skyrme to explain hadrons in nuclear physics \ skyrmions have been attracting consid- 
erable academic and technical interest because skyrmion presents transport properties associated with its 
global smooth topological spin texture^"^. Magnetic skyrmion and skyrmion crystal have recently been 
observed experimentally in chiral magnets^"^^, and have also been predicted to exist in coupled ordinary nano- 
magnets^^'^^. Furthermore, promising potential applications are envisioned for skyrmions in next-generation 
high-density storage and future microwave devices^^'^^ because of their small size of 10-100 nm and special 
topological nature. Considering the practical applications, the search for methods of manipulating skyrmions that 
are embedded in magnetic conductors or insulators (such as Cu20Se03) and the investigation of skyrmion 
dynamics are issues of vital importance. Spin-polarized current has been extensively used to drive skyrmions 
in conductors''^, but manipulation of skyrmions with a microwave field or a thermal current is still an important 
challenge^^'^^. Skyrmion dynamics under a microwave field presents various resonant spin-wave modes with 
small amplitudes and an alluring nonreciprocal transport effect^'^°. However, resonant excitation behavior of 
skyrmions in a microwave field has seldom been reported and investigating the nonlinear resonant excitation is 
helpful for understanding the dynamics of the skyrmion. Unlike a vortex with local deformation at the core and a 
bubble with local deformation at the domain walP^"^', a skyrmion moves with large non-local deformation. Based 
on its non-local deformation, the mechanisms of its dynamics are far from clear. Questions are what is the 
relationship between variation of deformation and trajectory of excitation, how can one control a skyrmion by 
resonance, etc. 

Here we report the first observation of flower-like resonant trajectories in Co/Ru/Co nanodisks under a single- 
frequency microwave magnetic field by means of micromagnetic simulations. The origin of the observed non- 
linear dynamics is interpreted as a combination of the large non-local deformation of the topological density 
distribution arising from the global smooth spin texture and the coupling of the resonant modes in coupled 
skyrmions. We find that the resonant excitations are controllable by tuning the frequency and the amplitude of 
the microwave field to adjust the deformation of the topological density distribution. The present work unveils the 
importance of the non-local nature of the topological density distribution for the skymrion dynamics and the 
prospects of skyrmions for spintronic applications. 
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Results 

A skyrmion is a nontrivial topological object with stable spin swirl 
and integer skyrmion number (topological charge) defined by the 
following formula^^'^^: 

Q= — qdxdy, q = m-(dxmxdym), (1) 

471 JJ 

where q is the topological density and m the unit vector of the local 
magnetization. Coupled skyrmions with opposite chiralities have 
been obtained in a Co/Ru/Co nanodisk without taking into account 
the Dzyaloshinsky-Moriya interaction (DMI), as shown in Fig. la. 
The topological charge of every skyrmion is minus one due to the 
same polarity. Figure lb presents the topological density distribution 
of a static skyrmion in the top nanolayer with a radial symmetry 
distribution. This global distribution is different from those of other 




Figure 1 | Sketch of the topological density distribution and two 
excitation modes of skyrmion dynamics, (a) Magnetization distribution 
in the top and the bottom nanolayers. The colors represent the magnitude 
of the out-of-plane magnetization component (m^) and the arrows 
indicate the direction of the in-plane magnetization component at every 
point, (b) The topological density distribution of a static skyrmion in the 
top nanolayer of a Co/Ru/Co nanodisk. q is the topological density, 
(c) Frequency dependence of the maximal radius R at which the guiding 
center can reach when skyrmion is excited by a microwave magnetic field. 
The field amplitude is set to be 100 Oe, except for the three points 
indicated by blue stars where 50 Oe is used because of the appearance of a 
very strong resonance. The excitation mode is clock- wise for field 
frequencies below 2.00 GHz (in the yellow region), and counterclock-wise 
above 2.00 GHz. In the green region, two modes with opposite rotation 
directions coexist. 



topological objects (such as Vortex, Meron and Bubble shown in 
Supplementary Fig. SI). In the center of the skyrmion, \q\ has the 
highest value. From the center to periphery, \q\ firstly decreases from 
8.8 X 10"^ to 3.2 X 10"^ nm"^, then increases to a local maximum 
near the skyrmion radius and finally gradually decreases to 0.4 X 
10"^ nm"^ at the boundary (\q\ as function of the radius is given in 
Supplementary Fig. S2). 

We investigate the resonant excitation behavior of coupled sky- 
rmions by applying a microwave magnetic field along the +x dir- 
ection on the top nanolayer, as shown in Fig. lb. If one of the coupled 
skyrmions is moved from its central equilibrium position by the 
torque of the external field, the other one is stimulated by the strong 
magnetostatic coupling and a non-local deformation of the topo- 
logical density distribution is created^'* (as will be shown later). 
Consequently, the trajectory of the skyrmion core becomes compli- 
cated and sometimes it is tricky to identify the exact core position. 
We have used the guiding center (Rx, Ry) of the skyrmion (the center 
of topological charge) to determine the skyrmion center position and 
the radius of the skyrmion motion. The guiding center is defined by^^: 



^ ^[xqdxdy ^ ^ $ yqdxdy 
^ $ qdxdy ' ^ ^ qdxdy ' 



(2) 



where q is the topological density defined in Eq. (1). Figure Ic shows 
the maximum steady-orbit radius R, which the guiding center can 
reach, as a function of the microwave frequency. The amplitude of 
the microwave magnetic field is set to be 100 Oe, except for the three 
points marked by blue stars where 50 Oe is used because of the 
appearance of a very strong resonance. In the microwave field, the 
guiding center of the skyrmion rotates around the nanodisk center 
and two resonance excitation modes are observed: mode 1 for clock- 
wise (CW) rotation of guiding center in the low-frequency region 
and mode 2 for counterclockwise (CCW) rotation in the high- 
frequency region, both of which are coupled with each other^^. The 
two modes with opposite rotation direction are analytic solutions of 
the generalized Thiele's equation. The forms of the solutions are^^: 



= -G/2M±J {G/lMy +K/M, 



(3) 



where G is the gyrocoupling vector, M the effective mass and K the 
stiffness coefficient. Equation (3) indicates that the high-frequency 
mode 2 is more sensitive to mass change than the low-frequency 
mode 1, which results in a much wider frequency range for mode 2 
than for mode 1, so that the frequency for mode fiipping is close to 
the low- frequency mode (around 2 GHz). Furthermore, in the 
Supplementary Fig. S6 the dependence of R on the field amplitude 
is also given. The frequency of the field is chosen to be 1.00 GHz, 
which activates linear dynamics: the plot of R as a function of the field 
amplitude shows an almost linear relationship. 

From the gyrotropic motion of the guiding center, we can obtain 
the eigenfrequencies of the coupled skyrmion, as demonstrated in 
Ref. 14. First, a square-wave pulsed magnetic field with 10-ns width 
and 50-mT magnitude is applied to the top nanolayer to shift the 
guiding center over some distance and then, after the field is switched 
off, the guiding center will gyrate around the center of the nanodisk. 
Eigenfrequencies of near 1.15 GHz (for a circular trajectory), and 
1.00 and 5.00 GHz (for a hexagonal trajectory) are obtained by 
means of the fast Fourier transform (FFT) spectrum of the trajectory 
of gyration. Figure Ic shows two pronounced peaks of R: one at 
1.13 GHz in the region of mode 1 and the other at 4.75 GHz in 
the region of mode 2. The asymmetry of the first peak may have 
by two origins: 1) the field amplitude used for the three star points is 
smaller than that for the other points; 2) the superposition of the two 
peaks with eigenfrequencies 1.00 and 1.15 GHz. The frequency of 
1.13 GHz of the first peak corresponds to the eigenfrequency of the 
circular trajectory. The frequency range of 4.60 ~ 4.85 GHz of the 
second peak, where a complicated dynamical behavior is found as 
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shown in Fig. 2, seems to deviate from the eigenfrequency of 
5.00 GHz for the hexagonal trajectory. Actually, the resonance peak 
at 5.00 GHz still exists, which indirectly may be concluded from the 
decreasing amplitude with increasing frequency which starts from 
5.00 GHz. This peak is more or less hidden behind a pronounced 
peak around 4.60 GHz. This frequency is an overtone of 1.15 GHz 
and activates the base mode, thereby giving rise to two excitation 
modes and a large amplitude. In the frequency range 4.60 ~ 
4.85 GHz, the rotation direction of the guiding center is initially 
CCW and becomes CW for steady orbits. The physical mechanism 
behind this novel non -linear dynamics and the appearance of the 
second peak needs further investigation. 

Now, we focus on the interesting frequency range 4.50 ~ 
4.85 GHz. The corresponding trajectories of coupled skyrmions 
due to resonant excitations are shown in Fig. 2. In a microwave field, 
the guiding center firstly rotates around the nanodisk center with an 
increasing radius until it reaches its steady orbits with periodic beha- 
vior (Fig. 2a). On the other hand, for quasi-periodic behavior, tra- 
jectories with the same radius have phase differences between 
different periods (Fig. 2b-2d). To illustrate the different results 
clearly, only the steady parts of trajectories are shown in the range 
of 30 ~ 50 ns. A steady orbit in a field with frequency of 4.50 GHz is 
circular with a radius of 12 nm, which is commonly observed in the 
dynamics of vortices and skyrmions excited by an alternating mag- 
netic field. However, when the field frequency is between 4.60 GHz 
and 4.85 GHz, the steady orbit transforms to a flower-like orbit with 
increased amplitude, which is quite unusual. The one-period traject- 
ories presented in the insets at the lower-right corners of Figs. 2b-d 
show five "petals". These petals evolve with time into ribbon- or 
flower-like trajectories with quasi-periodicity due to a phase differ- 
ence between the different periods. In addition, the phase shift 
increases from 0.02 to 0.15 radian if the frequency of the magnetic 
field changes from 4.60 to 4.85 GHz. The flower-like dynamics may 
be observed indirectly in an experiment: By using Lorentz transmis- 
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Figure 2 | Topological excitations of a skyrmion by microwave fields with 
different frequencies. Trajectories of the guiding center (R^, Ry) of a 
skyrmion in the top nanolayer in 30 ~ 50 ns for frequencies of the 
microwave field of 4.50 (a), 4.60 (b), 4.75 (c) and 4.85 (d) GHz. The 
amplitude of the field is 100 Oe. The insets in the upper- and lower-right 
corners show the corresponding FFT spectra and the corresponding 
trajectories in about one period, respectively. 



sion electron microscopy (LTEM)^°'^^ to obtain the magnetization 
distribution at different times and then calculating the position of the 
guiding center by means of Eq. (2). The FFT spectra of the traject- 
ories are shown in the insets at the upper-right corners of Fig. 2a-d. 
During the steady-motion time (30-50 ns), the FFT spectra do not 
change with the length of time. Only one resonant frequency is 
observed in the circular trajectory excited by a field of frequency 
4.50 GHz, whereas in a field of other frequency, an additional fre- 
quency around 1.15 GHz, which is reminiscent of the eigenfre- 
quency of 1.15 GHz, shows up. The FFT intensity of the intrinsic 
mode and the intensity ratio between the two modes change with 
frequency. 

To understand the flower-like trajectories better, a microwave 
field is applied to activate the topological excitation of the coupled 
skyrmion with a fixed frequency of 4.60 GHz but with amplitude 
varying from 20 to 500 Oe. The results are presented in Fig. 3. 
Surprisingly, when field amplitude is 20 Oe, the trajectory with a 
small radius of less than 2 nm is a circle instead of flower-like 
(Fig. 3a). The corresponding FFT spectrum shows only one peak at 
4.60 GHz. As the amplitude increases, the maximum radius 
increases, and trajectory changes to flower-like. With further increas- 
ing amplitude from 60 Oe to 140 Oe, the flower-like trajectory 
changes from a sunflower to a bauhinia, due to a reduced phase shift 
from 0.18 to 0.01 radian. As field amplitude reaches 500 Oe, the 
trajectory becomes quite distinct from the previous ones, with six 
"petals", instead of five, in each period (Fig. 3d). The corresponding 
FFT spectra in Figs. 3b-d are shown in the insets. Besides the peak at 
4.60 GHz, associated with the external field, another peak near 
1.15 GHz expectedly appears again. With increasing field amplitude, 
the intensity of peak at this frequency first increases and reaches its 
maximum at 100 Oe. Then the intensity decreases again with further 
increasing field amplitude, resulting in a decreased phase shift and 
different flower-like trajectories. It should be noted that, when the 
field amplitude is as large as 500 Oe, many frequencies apart from 
1.15 and 4.60 GHz appear, which accounts for the six "petals" of the 
trajectory. 
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Figure 3 | Influence of the amplitude of a microwave field on the 
trajectory of the guiding center. The frequency of the microwave field is 
4.60 GHz. (a), (b), (c) and (d) show the trajectories of the guiding center in 
30 ~ 50 ns when the skyrmion is excited by a field of 20, 60, 140 and 
500 Oe, respectively. The insets in the upper-right corners show the 
corresponding FFT spectra. 
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Figures 2 and 3 demonstrate that the flower-Uke trajectories are 
associated with an excitation of an intrinsic mode near 1.15 GHz, 
and the variation of the FFT intensity with the field ampUtude at this 
eigenfrequency leads to different trajectories. To understand why the 
1.15 GHz mode can be excited and exist stably in a microwave field 
of 4.60 ~ 4.85 GHz, Fig. 4 shows the topological density distribution 
of skyrmions corresponding to a maximum deformation in Fig. 3. 
The topological distribution of a static skyrmion has radial symmetry 
which, however, is broken as the skyrmion is activated by a micro- 
wave magnetic field. The deformation of the topological density 
distribution increases with field amplitude, becomes global if the 
amplitude larger than 20 Oe, and the topological density distribution 
becomes very complex when the field is as large as 450 Oe. To 
quantify the level of deformation, we start from Eq. (1) to calculate 
time-dependent value of the topological density change in a field as 
follows: 



D= J \q(t)-q(0)\dxdy. 



(4) 



Figure 5a shows the time dependence of D in the dynamical process 
of a skyrmion in a microwave field (4.60 GHz and 100 Oe). D varies 
periodically between 6. 1 and 2.6, in a similar manner as the trajectory 
of the guiding center in Fig. 2b, which implies a strong correlation 
between the topological density distribution and the non-linear 
dynamics of the skyrmion. Similar to the trajectory of the guiding 
center, also D is depending on the field amplitude. Figure 5b shows 
the field dependence of the average D (filled dots) and the oscillation 
of D (vertical bars) in one period during steady excitation. D is less 
than 0.5 at 20 Oe and initially increases monotonously to 4.50 at 
100 Oe. As the amplitude further increases, D starts to decrease till 
200 Oe and then increases again almost linearly to the highest value 
of 10.1. 

To unveil the role played by the energy terms behind the flower- 
like trajectory, we also show in Figs. 5c and 5d the energy variations 
corresponding to the D variation in Figs. 5a and 5b. To note here that 
all energy terms vary periodically with time but only energies that 
contribute dominantly to the total energy are shown in Figs. 5c and 
5d. The uniaxial anisotropy energy (E^n) and the demagnetization 
energy (Ejem) vary oppositely to each other with time and both of 
them are found to play a major role since the variation range of both 



the Ean and Edem is one order of magnitude as large as that of the total 
energy (Etot)- More importantly is the field dependence of the energy 
terms, which is similar to that of D. Figure 5d shows the field depend- 
ence of the average Etot, E^n, Edem (filled dots) and the Etot oscillation 
(vertical bars) in one period during steady excitation. The Etot oscil- 
lation, which shows the competition between sub-energy terms, 
exhibits a very similar behavior to the oscillation of D, suggesting a 
connection between D and competition among sub-energy terms. 
Moreover, the curves of Etot and Ean agree well with that of D but with 
Ean varying faster than Etot, whereas Edem has an opposite variation, 
all of which indicate that the competition existing between the Ean 
and Edem leads to the variation of D and, thus, to the flower-like 
dynamics. 

To clarify some possible issues such as whether the dynamics of a 
skyrmion in a microwave field can be simply depicted by the sky- 
rmion core displacement (m^ = — 1), we show in Fig. 6 the deforma- 
tion of magnetization distribution which corresponds to the 
topological density distribution of skyrmion in Fig. 4. From 
Fig. 6a, not only the core displacement (m^ = —1) but also the 
boundary deformation (m^ = 0) is observed in the magnetization 
distribution. As the field amplitude increases, both the core displace- 
ment and the boundary deformation become larger. To quantify 
such deformations, we show in Fig. 6b the maximum position dis- 
placement from the zero-field position of m^ (equals to — 1 and 0) as 
a function of the field amplitude. The displacement of core at m^ = 
— 1 increases initially to 18 nm and drops to 2 nm, and then rises 
again to 20 nm, which has a variation similar to D, whereas the 
displacement of the boundary at m^ = 0 increases all the way from 
0 to 17 nm as the field amplitude increases. It is interesting to see 
that, although the displacement of skyrmion core as a function of 
field amplitude is similar to that of the topological density distri- 
bution, the boundary deformation also has to be taken into account. 
Obviously, the dynamics of a skyrmion in a microwave field is more 
complex than that of a vortex and cannot simply be depicted by the 
motion of the skyrmion core (m^ = — 1). 

Figures 3-5 imply that it is the deformation of the topological 
density distribution that gives rise to the stable existence of the mode 
near 1.15 GHz and thus leads to the flower-like trajectories. 
Nevertheless, it is still unclear how deformation of the distribution 
leads to the stable existence of mode 1. To elucidate the role of the 
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Figure 4 | The topological density distribution of skyrmion in different microwave fields. The frequency of the microwave field is 4.60 GHz and the 
amplitudes are in the range 0 ~ 450 Oe. The figure shows the distribution at the time that D is maximal in one period. To compare the results easier, 
the distribution has been rotated by hand. The topological density distribution has radial symmetry without applied field and with increasing field 
amplitude, the deformation of the distribution increases and large non-local deformation takes place. 
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deformation of the topological density distribution, it is demon- 
strated in Fig. 7 that the flower-like trajectory can be tuned by con- 
trolling the level of deformation even at a frequency far from the two 
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eigenfrequencies (1.15 and 5.00 GHz)^^. The frequency of 3.45 GHz 
is chosen. The trajectory is approximately circular for the small field 
amplitude 100 Oe and, as shown in Fig. 7a, the corresponding FFT 
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Figure 6 | The magnetization distribution of a skyrmion in a microwave field. The microwave field applied has the same frequency (4.60 GHz) and 
amplitudes (in the range of 0 ~ 450 Oe) as that in Fig. 4. (a) The corresponding magnetization distribution to the topological density distribution of 
skyrmion (Fig. 4) under different field amplitude. The colors represent the magnitude of the out-of-plane magnetization component (m^) and the 
arrows represent the direction of magnetization projected in the plane at each point, (b) The maximum position displacement r^ax from the zero-field 
position as a function of field amplitude for both skymion core (mz = — 1 ) and boundary (m^ = 0). r^ax is chosen as follows. The initial position of m^ at a 
zero magnetic field is ro but changes to ri in a magnetic field, the maximum value of Iri — rol is selected for plotting. 
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spectrum has only one peak which is associated with the driving field. 
Changing the amplitude from 100 to 200 Oe to intensify the 
deformation of the topological density distribution leads to the 
appearance of a flower-like trajectory with an additional FFT peak 
at 1.10 GHz (Fig. 7b). A larger value of the amplitude seems neces- 
sary at 3.45 GHz than at 4.60 GHz to trigger the flower-like traject- 
ory of a skyrmion. These results further confirm that the deformation 
of the topological density distribution is vital for the stability of mode 
1 near 1.15 GHz and hence for the flower-like trajectories. 
Complementary to the trajectory of the guiding center in Fig. 7a is 
the trajectory of the average magnetization. Figures 7c and 7d show 
the corresponding trajectories of the average magnetization of 
coupled skyrmions in a microwave field with a frequency of 
3.45 GHz and amplitudes of 100 and 200 Oe. The three-dimensional 
(3D) plot of the magnetization shows periodic behavior when the 
field amplitude is 100 Oe. Its 2D projection on the xy plane is a circle 
while, when the field amplitude increases to 200 Oe, the 3D plot is 
complex with knots and shows quasiperiod behavior. Surprisingly, 
its 2D projection on the xy plane shows a flower-like trajectory, 
similar to that of the guiding center (Fig. 7d), which indicates that 
flower-like dynamics may be observed by detecting the volume- 
averaged magnetization signal at different times in the experiment. 

Figure 8 demonstrates the trajectory of the guiding center when a 
microwave field is applied to both the top and the bottom nanolayers. 
The field amplitude is 200 Oe and the frequency 4.6 GHz. Because of 
the opposite chiralities of the two skyrmions, the two guiding centers 
move oppositely and separate from each other. If there were no 
magnetostatic interaction between the two guiding centers, one 
would expect two independent microwave-field-induced flower-like 
trajectories in both the top and the bottom nanolayer. However, this 
center separation will change the interlayer magnetostatic inter- 
action which has an effect opposite to the microwave field and 
may destroy the flower-like trajectories. This is different from the 
case previously discussed with the microwave field only applied to 
one nanolayer, in which the magnetostatic interaction plays a pos- 
itive role and forces the other nanolayer to a synchronic flower-like 
trajectory. Though the two guiding centers move oppositely, the 
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Figure 7 | Influence of the topological density distribution deformation 
on the dynamics of a skyrmion. The trajectories of the guiding center in 30 
~ 50 ns are shown. The frequency of the microwave field is 3.45 GHz 
and the amplitude 100 Oe (a) and 200 Oe (b). The insets in the right 
corners are the corresponding FFT spectra. The corresponding trajectories 
of the average magnetization of the system in 100 Oe (c) and 
200 Oe (d) fields are also given. 
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initial rotation direction of the two guiding centers is the same 
(CCW), as shown in the upper-right insets. In addition, the magne- 
tostatic interaction between the two skyrmions is strong, which leads 
to a smaller distance of motion than that in Fig. 3. It should be noted 
that the rotation direction of both guiding centers does not change 
direction (indicated in the lower-right insets), which is different from 
the results in Figs. 2b and 3. 

Discussion 

We have observed complex flower-like dynamics with dual resonant 
modes activated by a single-frequency microwave magnetic field in 
coupled skyrmions by means of micromagnetic simulations. The 
occurrence of flower-like dynamics is associated with the large 
non-local deformation of the topological density distribution of 
coupled skyrmions. The novel non-linear resonant excitations can 
be controlled systematically by the amplitude and frequency of the 
external magnetic field, for example from a circular trajectory to 
different fiower-like ones. Because of its global smooth topological 
spin texture, a skyrmion can undergo a large non-local deformation 
of the magnetization distribution by a microwave field or by spin- 
transfer torque which leads to a global deformation of the topological 
density distribution. Therefore, we expect that similar fiower-like 
dynamics can also be activated in other skyrmion systems, such as 
a single skyrmion or a skyrmion crystal stabilized by the DMP°'^^'^^, 
not limited to coupled skyrmions, which may lead to observing 
fiower-like dynamics easier in an experiment. These findings con- 
tribute to an understanding of complex non-linear dynamics of topo- 
logical objects and are essential for manipulation of skyrmions in 
future applications. 

Methods 

Micromagnetics simulations. Three-dimensional (3D) micromagnetic simulations 
have been utilized to study the dynamical behavior of coupled skyrmions in Co/Ru/ 
Co nanodisks in an external magnetic field by means of the object oriented 
micromagnetic framework (OOMMF) code^^. For the cobalt material parameters, we 
used the saturation magnetization Ms = 1.4 X 10^ A/m, the exchange stiffness Agx = 
3.0 X 10"^^ J/m, and the uniaxial anisotropy constant K^ = 5.2 X 10^ J/m^ in the 
direction perpendicular to the nanodisk plane. The diameter of the Co/Ru/Co 
nanodisk was chosen to be 200 nm and the thickness of Co and Ru 18 nm and 2 nm, 
respectively. In the geometry, coupled skyrmions can be formed as a ground state^''. 
The cell size was 2X2X2 nm^ and the dimensionless damping a was 0.02. The 
interfacial coupling coefficient of the adjacent surfaces was —5 X 10"^ J/m^. The 
method of obtaining static coupled skyrmions is the same as that used in Ref. 14: 
Different initial magnetic states (vortex-like, in-plane-like, and out-of-plane-like 
initial states) were used to obtain the most stable state. To realize coupled static 
skyrmions in and experiment, we suggest to saturate the magnetization of the system 
by applying a magnetic field along +z direction (perpendicular to the disk plane), 
then reduce the field to zero and, finally, coupled skyrmions can be expected to form. 
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Figure 8 | Trajectories of the guiding center in a magnetic field with 
amplitude 200 Oe and frequency 4.6 GHz. The magnetic field is applied 
on both the top and bottom nanolayers. Trajectories of the guiding center 
in the (a) bottom and (b) top nanolayers. The insets in the upper- and 
lower-right corners are the initial and steady trajectories, respectively. The 
starting points are labelled by dots and the arrows indicate the rotation 
direction of the guiding center. 
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To investigate the resonant excitations of the coupled skyrmion, an external 
microwave magnetic field was applied along the +x axis to the top static skyrmion 
with different frequencies and amplitudes. So far, microwave technologies in 
nanostructures have achieved great progress, such as realizing spin waves and 
microwaves propagating in nanostructures via a nanostructured waveguide^""^^. In an 
experiment, we expect that the microwave magnetic field can be precisely applied to 
the top nanolayer by means of a nanoscale waveguide in future. 
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